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The time and frequency selective fading properties of radio channels may be characterized 
by evaluating the cross-correlation function between two received carriers as a function of 
their frequency separation. In practice such correlation functions have been measured 
using only the envelopes of the received carriers. Recent studies of the effect of time and 
frequency selective fading on digital data transmission have shown that envelope correlation 
information is insufficient for an accurate evaluation of system performance and that the 
complex envelope correlation function is needed. This paper presents an experimental 
technique for the measurement of the complex time-frequency correlation function of radio 
channels that uses independent frequency standards at transmitter and receiver. An ana- 
lytical study is made of the theory of operation of this system and of the effects of instabilities 
of the frequency standards. 



The transmission of a sinusoid through a fading 
radio channel results in the reception of a narrow- 
band random process. Let x(t) denote the received 
narrow band waveform resulting from the trans- 
mission of a sinusoid / c/s away from a "carrier" 
frequency of/ c/s. Then x(f) may be expressed as 



x(t)=B,e {T(f,t)e j2 " (f ° +f)t 



(1) 



where T(J, t), the complex envelope of the received 
narrow band process, is the Time Variant Transfer 
Function of the medium when measured with 
respect to the carrier frequency f . 

The correlation properties of this transfer function 
are defined by the correlation function 



R f .t{T,Q)=T*U,t)T(f+Q,t+T) 



(2) 



where * denotes the complex conjugate and the over- 
line indicates an ensemble average. Knowledge of 
the channel correlation function Rf,t(r, 12) is essential to 
determining the degradation of communication 
system performance caused by the time and frequen- 
cy selective fading of the channel. 

As defined in (2), the channel correlation function 
involves an ensemble average. Unfortunately, an 
ensemble of channels is not available and thus the 
ensemble averages may not be evaluated experimen- 
tally as averages over an ensemble of channels. In 
practice, the ensemble of interest is generated by 
considering the behavior of a single channel at a set 
of mutually exclusive time intervals of the same 
duration. In this way the behavior in a given time 
interval provides the typical member function of the 
ensemble, and ensemble averages are interpreted as 
appropriate time averages. 



The measurement problem is complicated by the 
fact that physical channels do not exhibit statistical 
regularity over indefinitely long intervals of time. 
Thus there is an inherent conflict between the desire 
to have a long measurement time in order to reduce 
measurement errors due to finite sample size and 
the desire to have a short enough measurement 
time so that the desired statistical parameters have 
changed negligibly due to the instability of the 
statistical structure of the channel. 

Bello [1963] has derived the quasi-WSSUS channel 
to account for the fading characteristics observed on 
many radio channels. A WSSUS (Wide Sense 
Stationary Uncorrelated Scattering) channel is one 
which, for analytical purposes, may be represented 
as a continuum of uncorrelated scatterers fluctuating 
with wide sense stationary statistics. A necessary 
and sufficient condition for a channel to be WSSUS 
is that the channel correlation function R f , t (r, 0) 
be independent of /, t; i.e., independent of the 
location of the time interval and the frequency 
band within which the measurement of the correla- 
tion function was made. Thus 



R f , i(t, Q)=B(t, 0) 



(3) 



for the WSSUS channel. 

The quasi-WSSUS or QWSSUS channel is one in 
which R f , t (r, 0) changes sufficiently "slowly" with 
changes in / and t. In particular, R f , t (r, 12) is 
assumed to change negligibly for values of f, t in 
a time-frequency rectangle whose time dimension is 
equal to the input signaling element duration plus 
the multipath spread of the channel and whose 
frequency dimension is equal the bandwidth of the 
signaling element. 



1 The work reported here was performed under subcontract to ITT Communications Systems, Inc., Paramus, N.J. 
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The utility of the QWSSUS channel concept is 
that, when applicable, it allows one to represent the 
channel as a WSSUS channel, at least as far as its 
effect upon any particular signaling element is 
concerned. 

Most radio channels exhibit a "fast" fading (e.g., 
a few c/s) superimposed upon a "slow" fading (e.g., 
hourly or diurnal variations). In order to measure 
the channel correlation functions by time averages 
the measurement time T m must satisfy the inequality 



T m « 



1 



e„ 



(4) 



where max is the bandwidth of the slow fluctuations. 
On the other hand, to keep the measurement errors 
due to finite sample size small one must have 



T m » 



B f 



(5) 



where B s is the bandwidth of the fast fluctuations. 
It follows that the short-time correlation functions 
can be measured only if the bandwidth of the fast 
fluctuations is a few orders of magnitude greater than 
the slow fluctuations. Fortunately, it appears that 
many radio channels exhibit this property. The 
following discussion assumes that (4) and (5) are 
satisfied and that over the measurement interval we 
can represent the channel by a hypothetical WSSUS 
channel. 

One further restriction needs to be imposed on 
the channel if the measurement techniques to be 
described are to succeed. This restriction is that 
the "spread factor" of the channel be less than unity, 
where the spread factor, k, is given by the ratio 



k B c 



(6) 



in which B c is the correlation bandwidth of the 
channel being measured. The reason for this restric- 
tion will become evident during the subsequent 
discussion. Fortunately, except for certain esoteric 
radio channels (e.g., the orbital dipole channel), the 
spread factor is less than unity and thus the measure- 
ment technique to be described is quite generally 
useful. 

We consider here the direct measurement of the 
correlation function R f >, ^(12, r), or rather R , t > (ft, r), 
since it is assumed that the center frequency j is 
chosen coincident with the "center" frequency of 
signaling waveforms to be transmitted. The basic 
measurement technique involves the transmission 
of two sinusoids separated by a frequency of ft c/s 
and the determination of the cross-correlation func- 
tion between the complex envelopes of the two 
received processes. Such correlation functions have 
never been measured for radio channels. The near- 
est approximation to this type of measurement is 
that performed by Claxton and Anderson [1962] in a 
lunar reflection study. They generated the two 
sinusoids as sidebands on an amplitude modulated 



carrier and determined the cross-correlation function 
between the envelopes (the magnitudes of the com- 
plex envelopes) of the two received processes. Thus, 
in effect, what they measured is given by 



Ba. l .(0,r)s|r*(-Q/2,^-T/2)||2 T (0/2,< / + r/2)| (7) 
rather than 



i?o.^(0,r) = T*(-0/2,f-T/2)T(ft/2,^ + r/2). (8) 

While the correlation function (7) is considerably 
easier to measure, it does not contain information 
due to decorrelation in the phase of the transfer 
function. Moreover, if the channel is indeed gaus- 
sian and QWSSUS, (7) does not provide sufficient 
statistical information for a complete characteriza- 
tion of the channel. Thus, for example, the error 
probability and intermodulation distortion calcula- 
tions performed by Bello and Nelin [1964a, b] cannot 
be carried through unless (8) rather than (7) is 
available. 

The measurement of (8) is difficult, not only 
because it requires that two narrowband processes 
at different center frequencies be manipulated so 
as not to introduce relative frequency instabilities, 
but also because it requires at the transmitter the 
generation of two sinusoids of different frequency 
with sufficiently small relative frequency instability. 

In figure 1 we present the block diagram of a 
measurement method which appears capable of 
measuring B(tt, r) for a range of situations of prac- 
tical interest. The success of the method depends 
upon the use of sufficiently stable clocks in the trans- 
mitter and receiver. If clocks of the required sta- 
bility are not available, a common timing source 
must be provided. Two methods for achieving a 
common timing source are the use of WWVL trans- 
missions and the use of a landline to provide the 
receiver with transmitter clock signals. 

We shall now discuss in some detail the theory of 
operation of the measurement technique of figure 1. 
A frequency standard in the transmitter whose out- 
put is nominally at F c/s is used in combination with 
a frequency synthesizer to generate a group of very 
stable sinusoids separated by multiples of F c/s. 2 
This frequency group is translated to an appropriate 
RF frequency f by heterodyning operations and 
then transmitted. Two types of frequency insta- 
bilities are inescapably associated with the trans- 
mitted output, namely, differential instability due 
to instability of the frequency standard and common 
instability due to instabilities in the heterodyning 
oscillators. From a mathematical point of view, 
common instabilities in the transmitter and receiver 
heterodyning oscillators may be regarded as part of 
the channel since they are also present when the 
channel is used for communication purposes. Thus 
common instabilities are not harmful provided they 
are small compared to the instability introduced by 



2 In practice it may be desirable to generate this frequency group by frequency 
modulating an oscillator with the frequency standard output and a modulation 
index chosen to given sidebands and carriers of not too different amplitudes. 
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Figure 1. Block diagram of complex correlation fundi on measurement method. 



the medium itself — a situation readily achieved in 
practice. On the other hand, differential instabilities 
are quite harmful since they are falsely interpreted 
by the measurement equipment as a frequency se- 
lective behavior. Moreover, this frequency selective 
behavior caused by differential frequency instability, 
in contrast to the effects of common frequency insta- 
bility, are not present when the channel is used for 
communication . 

From the above observations of equipment fre- 
quency instability, we conclude that in our analysis 
of the measurement technique of figure 1, it is 
sufficient to neglect the common frequency instabili- 
ties introduced in the transmitter and receiver, 
assuming that they have been included as part of the 
transmission medium. Moreover, as may readily be 
appreciated, it is only the relative instability between 
the transmitter and receiver frequency standards 
that enters in the accuracy of the measurement 
technique, so that for analytical purposes it is suffi- 
cient to assume a perfect frequency standard at the 
receiver and a hypothetical frequency standard at the 
transmitter whose frequency instability is that ob- 
served for the relative frequency of two independent 
frequency standards. 

In view of the above, the transmitter output y(t) 
may be expressed as 



2/(0=S cos [2TrUo+nF)t+n5Ft+ncj)(t)] (9) 



where /'o is the RF frequency, 8 is a random variable 
describing the long term percentage drift between 
two independent frequency standards, and <£(£) 
describes the short term phase fluctuations about 
the long term phase drift 8Ft. An initial analysis 
was carried through, including the phase fluctuation 
term <t>(t). However, it was found to produce negli- 
gible measurement error compared to the long term 
phase drift SFt. Thus, to simplify our discussion 
we will neglect the short term phase fluctuations and 
use 

y(0 = S cos [2ir(J +nF)t + n8Ft + n6] (10) 

in place of (9), where 6 is a constant denoting a 
random initial phase at t=0. 

According to our definition of T(f, t) we see that 
if y(t) in (10) is the transmitted process; then, after 
translation to an IF frequency of fi c/s, the received 
process complex envelope is given by 

w(0 =S TinF, t)ef 2Hfl+nF+nF8U+ne . (11) 

To be exact one should use T(nF-\-nF8, t) in the 
summation in (11) rather than T(nF,t). However, 
for any physical channel these will be indistinguish- 
able since 8 will be extremely small (10~ 10 perhaps). 
Each of the N+l fading carriers at the frequencies 
fi+nF; n = 0, . . . , N is extracted by means of a 
filter bank. Some attention must be given to the 
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design of these filters so that, with respect to their 
center frequencies, their amplitude and phase char- 
acteristics are essentially identical over the band- 
width of the input fading carrier. This should not 
be a very difficult problem since for virtually all 
channels of interest the required value of F to sample 
the channel correlation function accurately will be 
much larger than the fading bandwidth, and thus 
each filter bandwidth can be much wider than the 
bandwidth of the process it is desired to filter. 

At the receiver a frequency standard and synthe- 
sizer generates a set of sinusoidal components at the 
nominal frequencies /, +nF; n=0, . . . N. As dis- 
cussed above, we have analytically taken account of 
receiver frequency standard instabilities by including 
them in the transmitter frequency standard insta- 
bility and heterodyning oscillator instabilities by 
including them in the channel. Thus the output of 
the receiver frequency synthesizer in figure 1 may be 
represented as the set of discrete frequency compo- 
nents cos 2ir(f 1 +nF)t; n = 0, . ..,2V. With the 
aid of 90° phase shifters, the set of components sin 
2ir(fi+nF)t; n = Q, . . ., N is also generated, and 
then coherent detection of the in-phase and quadra- 
ture components of the various received carriers is 
performed. The bandwidths of the low-pass filters 
following the multipliers are just wide enough to pass 
the fluctuations due to the fading of the correspond- 
ing carrier. It is readily seen that the outputs 
X n (t) and Y n (t) in figure 1 are given by 



t + nO] 



X n (t)=Re{T(nF,t)e il2 "' n 



F„(0=lm { T(nF, t)e n2rnFSl+ne] 



(12) 



It will now be demonstrated that by appropriate 
processing of the set of waveforms {X n (t), Y n (t)} it 
is possible to determine the required channel corre- 
lation function. The set of waveforms [X n (t) , Y n (t) } 
may be recorded on tape or film and later processed 
by a digital computer or else processed in real time 
with an analog computer. 

As discussed previously, we shall assume that over 
the measurement period we may assume the channel 
to be WSSUS. Thus the correlation function 



T*(f,t)T(J+Q,t+r) 

<T*(f, t)T(f+Q, t+r)>=B(Q, r). (13) 

If we form the functions X 0n (f), Y 0n (t) given by 

XoM^Xo^X^t + rO-YoiDY^t + rO 

Y 0n (t)=X (t)Y n (t + r 1 )+Y () (t)X n (t + r l ) (14) 
we readily determine from (24) that 
X 0n (t) = ~Re { T*(0, t)T{nF, t+r l )e i ^ n9n+n ^ } 
Y 0n (t)=Im { T*(0, t)T(nF, t+ Tl )e" 2 * n8Ft+ ^ } (15) 



where in (15) differs from in (12) in containing 
an additive term 2im8FT l . By performing the short 
time averages 

I rt+T m i2 
{X Qn {t))T m =- Fi r I m , o X 0n (s)ds 

(Y 0n (t))T m =^r f™ Y 0n (s)ds (16) 

we can determine real and imaginary parts of 

p in6 rt+T m i2 
Z 0n (t)~ " T*(0,s)T(nF, s+rje^'ck 

=^ f Tm/2 r*(0, s-t)T(nF, s-t + Tl )e MF °ds 

■J-m J ~T m /2 

(17) 

where 6 in the second equation in (17) differs from 
the first equation by the quantity —27rn8FT m . 

Now the product T*(0,s—t)T(nF,s—t+r) contains 
a direct current value plus a fluctuating component ; 
i.e., 

T*(0, s-t)T(nF, 8-t+ri)=R(nF, n) +/,(*) (18) 

and it is the purpose of the time average to average 
out the fluctuations j n {t). Assuming that T m is 
large enough to perform this averaging we find that 

Z Q ^)=e^[R(nF,rO+Ut)}e^nT m 
~e jnd R(nF, n){e^ nm )T m 
= e ind sine [nF8T m ]R(nF, r x ) (19) 



where 



sin tx „ (irx) 2 . (ttxY 
smcx= =1- 



ttx 



(20) 



3! ' 5! 
From (19) we note that 

\Z 0n \ = \smcnFST m \\R(nF, n )\ 

^Z 0n =n9+^B(nF,n). (21) 

Thus the clock frequency instabilities appear as an 
error only in the measurement of the magnitude of 
the required channel correlation function while the 
ambiguity in initial phase affects only the angle of 
the channel correlation function. The greatest 
error in the magnitude occurs at n=N. It is 
probably sufficient for most purposes to require 



l-|sinciVF5rJ<0.1 
then, from the expansion in (20) 
NFST m <0.25. 



(22) 



(23) 
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To obtain an idea of the long term stability 
required in a rather extreme case assume 



NF=10 megacycles per second 
I 7 * = 100 seconds 



Then from (23) we find that 

<5<2.5X1(T 



(24) 



(25; 



which is just within the capability of crystal con- 
trolled frequency standards. 

We will now describe a method to remove the 
effects of the angle 6. Unfortunately this method 
not only removes the effect of the angle 6, but it 
removes the effect of any linear variation of <£i?(ft, r) 
with 12. This fact becomes obvious when we write 
nO in the form 

nd=n2irFT 1 (26) 

where 



Zi= 



2tF 



(27) 



It is readily demonstrated that a linear variation of 
<£i?(0, r) with Q of the form 27rl2r is directly inter- 
pretable as the existence of a mean path delay 
t . Thus, from (26) we see that the presence of the 
angle nd may be interpreted physically as adding to 
the channel an unknown path delay T x . 

For purposes of discussion, assume that any mean 
path delay has been absorbed by T x so that<£i?(12, r) 
contains no linear variation with 12. In actuality 
the phase of Zo n will be measured only in a range of 
2ir, say from — rr to it. Thus let 

<f>n={<Z 0n } Modulo 2tt (28) 

be the angle actually measured for -£Z 0n . Form the 
second order difference 

A 2 (/>,-{[(/> w -0«,_i]-[0n-i-^-2]} Modulo 2tt. (29) 

Then it is readily demonstrated that 

A 2 cj> n =A 2 [^R(nF, n)] Mod 2tt. (30) 

If the procedure is now reversed and it is assumed 
that <£i?(12, n) varies small enough in an Q interval 
of duration F so that 

A[<5(^r 1 )]=^«(^r 1 ) 

-^R(\n-l]F,r 1 )«27r (31) 

one finds that the resulting set of points will form 
a smooth curve except where it reaches the limits 
— 7r or 7r. These points of discontinuity may be 



removed by adding or subtracting 2x and continuing 
the curve smoothly. Presuming this has been 
done, the resulting curve will differ from the actual 
<£i?(Q, ri) only in a linear variation of the form 
27T127 7 ! and in a constant value which is equal to 
<£//(0, n). This initial angle may be found from 
</ o since 

<£Z o=<R(0,T 1 ). (32) 

As a final note it should be pointed out that if 
one may assume that the channel has gaussian 
statistics it is possible to determine \R(nF, n)| by 
an alternate procedure that is unaffected by the 
instability of the frequency standards. This method 
involves a determination of the set of correlations 
between envelopes of the outputs of filter bank in 
figure 1; i.e., a determination of 



W« n =(\T(S),t)\\T{nF,t + n)\)T m . 



(33) 



If one assumes that T(f, t) is normally distributed, 
then, assuming that T m is large enough to average 
out the fast channel fluctuations one may show from 
the properties of a gaussian process that 



w 0n 



#(0,0) 
where 



2E[ P (.nF, Tl )]-[l-p(nF, ^IKlpiflF, r,)] 

(34) 



p(nF, n ) = 



R(nF, Tl ) 



K(0, 0) 



(35) 



and E[ ], K[ ] are complete elliptic integrals of the 
second and first kind, respectively. Because of 
the monotonic relationship between W JR(0, 0) 
and p(nF, rO, it is possible to determine \R(nF, n)l 
once W Qn and R(0, 0) are measured. 
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